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SUMMARY 


This paper presents the results of a parametric study of nine 
analytic models that are commonly used for representing aerosol size 
distributions; discusses methods by which best-fit estimates of these 
parameters can be obtained; describes a catalog of graphical plots, 
depicting the parametric behavior of the functions; and explains pro- 
cedures for obtaining analytical representations of size distribution 
data by visual matching of the data with one of the plots. These 
analytic models consist of mathematical functions, with up to four 
adjustable parameters, which are referred to as follows: power law; 

regularized power law — in which singularities that occur in a simple 
power law as radius goes to zero are avoided; modified gamma distribu- 
tion; inverse modified gamma distribution — in which the variable is the 
inverse of radius; log normal distribution; normal distribution; 
generalized distribution function; and, power law generalized distribu- 
tion function. The mathematical properties of these distribution 
functions, such as, mode radius, modal value, limiting behavior of the 
distribution as radius goes to zero or infinity, moments of the dis- 
tribution, and, the dominant role played by each parameter in the 
model. These properties are summarized in a table in the report. The 
model catalog consists of log-log graphs (6x4 cycles) and/or in 
some cases, semi-log graphs on which the model size distributions are 
separately plotted. For each of the nine models, sets of plots have 
been produced, each set depicting the behavior of the model as a 
function of only one of its parameters, the other parameters being 
kept constant. Thus, each catalog graph contains several plots, one 



2 


corresponding to each value of the parameter being varied, the para- 
meter values being chosen to cover the range of variations in aerosol 
size distribution likely to be encountered. 

The catalog is meant to assist researchers in obtaining an 
analytic representation for their empirical size distribution data. 

This is done by means of a relatively quick visual fit to one of the 
plots in the catalog; and, if the fit is not good enough, which is more 
likely to be the case, then use the parameters corresponding to the 
closest plot as initial estimates in a least-squares type program to 
obtain the best fit to the data. For this purpose, the experimental 
data should be plotted on transparencies provided with the catalog. 

The closer the initial estimates are to the best-fit values, the 
better the assurance of reaching a fast convergence. 

In addition, it is shown that the same experimental data can often 
be represented with equal accuracy by more than one analytic function. 

INTRODUCTION 

Even though atmospheric aerosols are known to possess a variety of 
shapes, the description of their physical structure is immensely, simplified 
if they are assumed to be spherical. The size spectrum of atmospheric 

aerosols is, in general, continuous and covers over four decades in radii, 

_3 

viz, 10 to 20 pm (Ref. 1). Of the basically four ways in which the 
empirical size distribution (SD) data can be represented, namely, tabular, 
histogram, graphical and analytical (Ref. 2), the last one is usually 
employed due to the fact that there exist regularities in the gross 
structure of atmospheric aerosols which exhibit behaviour similar to that 
of a variety of mathematical functions. An analytic function generally 
encompasses in a smooth way the main features of the aerosol physical structure 
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While admittedly unrealistic in its smoothness, the analytical, represen- 
tation has the following advantages, namely those of convenient 
adjustability to obtain a best fit to the experimental data; compact 
representation of the dependent variable (the SD) in the form of 
estimated parameters of the fitted distribution; construction of 
reasonable and convenient models; and, carrying out systematic, para- 
metric modeling of the aerosol optical properties. 

The success of the analytic representation approach depends upon the 
the selection of an appropriate mathematical function to approximate 
the actual size distribution data. This may not always be possible 
(Ref. 2); often a linear sum of mathematical functions may provide a 
good representation. Thus, there seems to be no "special” analytic 
function that can be said to be unique in representing aerosol SD's. 

The choice of the function is to some extent dictated by the modeler’s 
taste. However, ultimately, it is only when the fitted analytic function 
leads to results that closely fit the experimental optical (scattering/ 
extinction) data and at the same time falls within the typical physical 
domain of atmospheric aerosols, that the analytic function may be assumed 
to represent the aerosol SD. 

The main purpose of this paper is to describe the use of a catalog, 
depicting graphically the parametric behavior of some analytic functions 
with up to 4 parameters, that could assist researchers in obtaining analyt- 
ical representation of their empirical particle size distribution data. 

The terminology used, the distribution functions considered and the basis 
for their selection, and a description of the various models and their 
mathematical properties precede the discussion on the catalog. 




4 


The valuable discussions in connection with this work and the careful 
review of the paper by M. A. Box, Institute for Atmospheric Optics and 
Remote Sensing, and W. P. Chu and K. H. Crumbly, NASA-Langley Research 
Center, are gratefully acknowledged. 

Terminology for Aerosol Size Distributions 
For clarification, the terminology for the particle SD’s used is dis- 
cussed in this section. The physical structure of aerosols (atmospheric or 
artificial) can be represented, in general, in terms of the number, area, 
volume or mass of aerosol particles per unit volume per unit radius at r. 
However, in this paper the discussion will be restricted to the SD represen- 
tation in terms of the particle number and radius, so that only the radius-: 
number distribution, log radius-number distribution and cumulative size dis- 
tribution will be considered. Their definitions, used in this report, follow 
those given in Ref. 3. (The dependence of the SD on factors such as altitude, 

composition, etc., will not be considered here.) 

-3-1 

The radius-number distribution n(r) , (cm ym ), is defined as the number 

3 

of particles per unit volume (cm ) within a unit radius range at r measured 
in ym. Thus, 

3 “1 

n(r) = -dN ( r ) / dr = dN^r) / dr, (cm“ ym ) 0-) 


where N(r) (N^r)) is the cumulative oversize (undersize) distribution 
function: (Ref. 4) • 

—3 

The log radius-number distribution n L (r) (cm ) is defined bv 
n L (r) =-dN(r) / d(log 1() r) - 2,3026 r n(r), (cm" 3 ) , 
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Junge found it convenient to handle the wide range of atmospheric aerosol 

size distribution data by plotting n^(r) as a function of r, on a log- 

log scale. This method of plotting has the advantage that it represents 

the particle concentration as well as the size distribution. 

The cumulative size distribution represents the total number of 
3 

particles per cm that have radii greater (less) than r, is represented 
by N(r) (N u (r)) and is called the cumulative oversize (undersize) distribu 
tion function (Ref. 4). In essence 


N(r) = 


n(r")dr" 


-3 

n (r') d log (r'), (cm ) 
log 10 (r) L 10 (3a) 


n (r) 
u 


J 


n (r " ) dr " = 


( log i0 (r) -3 

n (r*) d log (r'l, ( cm ) 

X (3b) 


In this paper, only the cumulative oversize distribution (COSD) will be 

-3 

discussed. Note that N(0) would then be the total number density (cm ) 
of the particles. 

Distribution Functions and Selection Criteria 
Given some empirical aerosol size distribution data, the problem 
is to find an analytic function that will most closely represent this 
data. Examples of mathematical functions of up to two parameters are 
the normal, gamma, binomial and exponential distribution functions; and 
of those with more than two parameters are the Weibull, Johnson and 

Pearson distribution families. These distributions are discussed in 

■ ! 

detail in many books on prcjbabi lity. (See Ref. 5.) These functions admit 
almost every type pf probability distribution, except composite distri- 
butions made up of several (distinct populations, such as 
multimodal distributions. In addition, there is another versatile 
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distribution, referred to as the generalized distribution function (GDF) 
(Refs. 6 and 7) which is derived from the Wood-Saxon function. 

In the selection of an analytic function to represent the size 
distribution n(r), the following criteria must be taken into account: 

1. The function is not singular for 0 ^ r ^ 

2. It is easily integrable over r; 

3. It can represent the main features of the gross structure of 
the aerosols by a minimum number of adjustable parameters. 


Symbols 


DN(R) /DLOGIOR 
DN(R)/DR 

\ 

n(r) 

n (r) 

j 

n L (r) 

N(0) 

N(r) 

N(>r) 


[N u (r)] 


P i' X 


P i 


= 1 , 2 , 3 , 4 ,, 


P i 


P 


24 


•42 


LIST OF SYMBOLS AND ACRONYMS 


adjustable constants (Defined in Eq. 18a) 
average particle surface area 

n (r) = dN(r)/log r; notation used in the plot labels 
L 10 

n(r) = d N(r) /dr; notation used in the plot labels 
moment 

r -3 -It 

radius-number size distribution Lem ym J 
n(r) for j th mode (Defined in Eq. 18a) 

-3 

log radius-number size distribution [cm ] 

, -3, 

number density Lem J 

cumulative oversize [undersize] distribution [cm ] 

same as N(r) ; notation used in the plot labels 

adjustable parameters in mathematical functions; 
given values of parameters 

estimated values of parameters 

reestimated values of parameters 

(P 2 + 1)/P 4 

(p 2 ‘ 1)/P 4 
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r 

r 

r . 

3 

r 

m 

r r 

Lm 



V 

Acronyms 

CG 

CP 

COSD 

DT 

GD 

IMGD 

LND 

MGD 

MODELR 

ND 

NLLS 

PLD 

PLGD 

QCM 

RM 

RPLD 

SP 

SD1 , SD2 


particle radius 

average particle radius 

j value of r 

mode radius for n(r) 

mode radius for n (r) 

L 

limits of r 

average particle volume 

catalog graph 
catalog plot 

cumulative oversize distribution, N(>r) 

data transparency 

generalized distribution 

inverse modified gamma distribution 

log normal distribution 

modified gamma distribution 

represents model number 

normal distribution 

nonlinear least squares 

power law distribution 

power law generalized distribution 

quartz crystal measurement 

r , the mode radius 
tn 

regularized power law distribution 
size distribution 
terms 1 and 2 of SD 
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SIZE DISTRIBUTION MODELS AND THEIR MATHEMATICAL PROPERTIES 
Analytic models suitable for representing aerosol size distributions 
include the following mathematical functions: 

1. Power Law Distribution (PLD) 

2. Regularized Power Law Distribution (RPLD) 


3. Modified Gamma Distribution (MGD) 

4. Inverse Modified Gamma Distribution (IMGD) 

5. Log-normal Distribution (LND) 

6. Normal Distribution (ND) 

7. Generalized Distribution (GD) 

8. Power Law Generalized Distribution (PLGD) 

The expressions and the mathematical properties of the functions 
will be described in this section. Here the model distributions represent 
the radius-number distribution n(r) , from which the corresponding expres- 
sions for n^(r) and N(r) are derived. The properties of interest are: 

the mode radii, for n(r) and n (r) ; lower limit, asymptotic, and para — 

L 

tin 

metric behavior of the functions; and the k moment of the models. 

The mode radius r for n(r) is given by the solution of 
m 


d n(r) 
dr 


= 0 


and that for n L (r), by the solution of 


d n L (r) 
dr 


= 0 


th 


The k moment is given by 


\ = 


r n(r)dr 


'0 
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The moments are useful for calculating properties of the distribution 
such as number (N(0)), average radius (r) , average area (A), and 
average volume (V) of aerosol particles in a unit volume, as shown 
here: 

N(0) a'M 
r = M 1 /M 0 

A = ttM 2 /M 0 

These properties for the eight size distributions models for n(r) and N(>r) 
are summarized in Appendix I for quick reference. 

A description of how one can determine good estimates of the para- 
meters from the gradients, modal values and moments of the models will be 
given in a later section. In all these models, the adjustable parameters 
are represented by p^, p 2 , Pg > •••♦, where p^ is the scaling parameter and 
is chosen so that the maximum value of the function is unity. 

Model 1 — The Power Law (PL) Model 
This model, known as the Junge power law, was proposed by 
represent his continental aerosol Sp data and is given by 

n(r) = p 1 r , r - r - r 2 
or, alternately, by 

( 1 ~ p 2 ) < < 

n L (r) =2.3 p^r , r^ - r - r 2 

It has a COSD of the form 


Junge to 


(4a ) 


(4b) 
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N(r) = — (r 


P 1 , 1_P 2 1_P 2v < < 

- r 2 ) , r i s r * r 2 


P 2 -! 


(4c) 


Junge used O.Olym for and l.Opm for r 2 but other values could be 


used . 


fctl 

The k moment for the distribution is given by 


l+k-p„ l+k-p„ 

\ = P l (r l “ r 2 ) / (P 2 ~ k-1) 


(4d) 


This model becomes singular at r = 0, if r 1 = 0; has all its moments 
infinite if r^ = 0 and r 2 = oo ; and has no mode radius ( r m )* Even though 
this model may not always represent a real situation, and does not meet 
the selection criteria, it is popularly used as it readily gives analytically 
tractable results. The model is graphically presented in Figs. (1A.1 - 1C.1). 

Model 2 — The Regularized Power Law (RPL) Model ■ 

In order to eliminate the singularity at r = 0 that occurs in Model 1, 
without losing its power law behaviour at large r, one may use a regularized 
form of the power law, 

(5a) 


The mode radius is given by 


n(r) = 


f P !l 

( r/B2 ) 

- P 2 > 

r, . . , 


P 3 _1 


r 

m 


= P, 


P 3 -i 


(1 + P 3 (P 4 -1)J 


(5b) 
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and the maximum value is 


1 - 


1 +P 3 (P 4 - 1 ) 


n(r ) = - 

m p 


1 3 ( 1 + P 3 (P 4 - 1 )) 

2 ( , P 4 

P 3 P 4 


(5c) 


The limiting behavior of the function is as follows! 


P 1 P 3 ~ 1 

As r Q, n(r)* — — (r/p„) 

p 2 


(5d) 


and, as n(r) ^ (r/p 0 ) 


P 1 -(1 + P 3 (P 4 -1)) 


(5e) 


The log radius number distribution is 


n L (r) = 


2.3 p 1 (r/p 2 ) . 

'p 3 p, 
(l+(r/p 2 ) J ) 4 


( 5 f ) 


Its maximum value occurs at 


r lm = p 2 (p 4 -1 } 


■1/p, 


(P 4 -l) 

so that n T (r. ) = 2.3 p, *— 

u im 1 p 


p 4 -1 


(P 4 ) 


(5g) 


The limiting behavior for this form of the distribution is as follows 


As 


r 0 , n(r) =; 2 .3 p^ (r/p 2 ) 


(5h) 


and, as r + 00 , n(r)~ 2.3 p x (r/p 2 ) 


-P 3 (P 4 -1) 


(5i) 
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The COSD is given by (Ref. 8, No. 3.194/2 and 9.121/1), 


N(r) = 


P 3 <P 4 ’- 1 ) P 3 P 4 ~ l 

J 4 (1 +(r/p 2 ) J ) 4 


(5j) 


The limiting behavior for this function is as follows, 


As r-?- 0, N(r) = - (1 - (p^-l) (r/p 2 > J ) 

-P 3 ( P 4 - 1) 


and, as r ->■ 00 > N(r) ~ — — -ppr (r/p„) 

P 3 vP^ - a / 4 


(5k) 


(5 A) 


The moments for the RPL are given by (Ref. 8, No. 3.194/3 and 
Ref. 9, p. 103) , s;; 


M, = 

k p 


p l p 2 ni+k/p 3 )r(p 4 -i-k/P 3 ) 


r (p 4 > 


, k <c p 3 (p 4 -l) 


(5m) 


where F is the complete gamma function. 

The parameter p 2 has the main effect on mode radius, being a multi- 
plicative factor, while p 3 and p^ control the positive and negative gradients, 
and hence polydispersity . The parameter p^ controls the positive gradient 
while both p^ and p^ influence the negative gradient. The model is 
graphically presented in Figs. (2A.1 - 2C.4). 


Model 3 Modified Gamma Distribution (MGD) Model 

Model 3 is the modified gamma distribution function. Deirmendj ian 
(Ref. 10) has shown that this function can be used to describe various 
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types of realistic aerosol distributions. For instance, by assigning 
different values to the parameters P 2 and p^ one obtains models such as 
Haze H, Haze M, Haze L, Cloud C3, etc. 

The radius~pumber distribution is given by 


P 2 

n(r) = p 1 r 


exp (-p 3 r P 41 


(6a) 


Its mocje radius is 


r 

m 


Pi 


P 3 P 4 


1/p, 


(6b) 


and its maximum value is 

f Po 


n(r) = p 
m I 


' p2 /p 4 


P 3 P 4 


exp (-P 2 /P 4 ) 


(6c) 


The limiting behavior of the distribution is as follows 

Po 


As r -*■ 0 , n(r) 


Pi r 


(6d) 


and, as r + <*>, n(r) ■* 0 as p L exp (-p r 4 ) 


(6e) 


The log radius-number distribution is given by 

p 2 +1 P 4 

n (r) * 2.3 PjT exp.(-p r ) 


(6f) 
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It has a maximum value at 



( P24 /p 3 ) 



where 


p 24 ~ ^ p 4 


(6g) 


so that 


24 


n L (r m ) = 2 * 3 Pi <P 2 4 /p 3 ) exp ( " p 24 ) 


(6h) 


The limiting behavior of the function is as follows, 


As r ->0, n L (r) ^ 2.3 r 


■V 1 


(61) 


and, as r -* oo, n^(r) 0 as 2.3 exp(-p^r ; 


The COSD is given by (Ref. 8, No, 3-381/3), 


(6 j ) 


N(r) * Pl r (p 24 , P 3 r *) / p 4 P 3 24 P 3 > 0 


(6k) 


The limiting behavior for the COSD is as follows (Ref. 11, 6.5.3, 
6.5.12, (j .5 .32) , 


As 


p 24 ' p 2 +1 


r -v 0, N(r) « p^r (p 24 ) - P 3 r /p 24 ) / P 4 P 3 


24 


(6JO 


p P 2 -P,+! P 4 

and, as r -»■ oo, N(r)~ r exp(-p r ) 


P 3 P 4 


(6m) 
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The moments for the distribution are given by (Ref. 8, No. 3.381/4), 

“(PoA+k/p.) 

\ " ^ p l^ p 4^ p 3 ^ p 24 +k/ ^ P 4^ ’ p 2 +k+1 > 0 

The parameter p^ controls the mode radius and the parameters p ^ and 
p 4 control the polydispersity . The parameter p 2 determines the limiting 
behavior as r+O while the parameter p^ determines the limiting behavior 
as r •+ 00 . The model is graphically presented in Figs. (3A.1 - 3C.5). 

Model 4 — Inverse Modified Gamma Distribution (IMGD) Model 

This distribution has the same form as Model 3 except that the inverse 
radius is used. This results in an exponential fall-off at the small size 
end and power law behaviour at the large-radius end. Twomey (Ref. 12) suggests 
this form of the modified gamma distribution for dry aerosols. 

The radius number distribution is given by 


P 4v , p 2 


n(r) = p x ex;p(-p 3 /r ) / r 


(7a) 


Ita mode radius is given by 


r P 3 P 4l 


1/p 


4 ■ 


m 


( 7b) 


and the maximum value is 


P. 


p 2 /p 4 


n(r ) 
m 


Pi 


P 3 P 4 


exp(-p^/p 4 ) 


( 70 
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The limiting behavior of the distribution is as follows, 

p 4 

As r ^ 0, n(r) 0 as expC-p^/r ) 

~P 2 ' : 

and, as r -»■ oo, n (r) ~ p r 

The log-radius number distribution is given by 

p, (p„-l) 

n fj (r) = 2.3 p^ exp(-p^/r ; / r 

Its mode radius is given by 

.■ 1/p. ■. ■ ■■■; 

r lm "" ^ P 3^ P 42^ ’ p 42 = ^ p 2 “ 1 ^ p 4 

and the maximum value is 

“L^liu 5 = 2,3 p l 42 exp(-p 4 . 2 ) 

The limiting behavior of the distribution is as follows, 

p 4 

As r ->0, n^(r) -► 0 as 2.3 p^ exp(-p^/r ) 

-(Po-D 

and, as r -* °o, n L (r) ~ 2.3 p t r 

The COSD is given by (Ref. 8, No. 3.381/1) 

K(r) ■ | * <I>42-»3 r ’ P4 >" *42 >0 ■ 

where y is the incomplete gamma function. 


(7d) 

(7e) 

( 7 f ) 

(7g) 

(7h) 
( 7 i) 

< 74 > 

(7k) 
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The limiting behavior of this distribution is as follows (Ref. 11 ; 
No. 6.1.1, 6 ,5 .£ , 6.5.12). 


-P 


As 


ri ? Pa P/(P A9 ~ ^ 

r + 0, N(r) =? -I [p 3 r (p^ 2 ^ - exp(-p 3 /r 4 ) / p 3 r 4 


] ( 7 £) 


and, as x 00 , N(r) ~ r 


-(p 2 -l) 


/ (p 2 -l)' 


(7m) 


The moments of the distribution are given by (Ref. 8, No, 3,381/4) 


pi -<? 42 -i>'/p 4 ) 

k * p 4 p 3 


r (p 42 -k/p 4 ), k < P 2 -1 


(7 n ) 


The parameters p 2 and P 4 control the rate pf fall-off at large end 
small radii, respectively, and hence control the polydispersity. The 
parameter p^ controls the mode radius. This model is graphically pre- 
sented in Figs. (4A.1 - 4C.4). 

Model 5 — The Log Normal Distribution (LND) Model 
The log normal distribution generally provides a better description of 
particle size distribution than the normal distribuiton (discussed later on) 
because particle sizes, like many naturally occurring populations, are often 
^symmetric, In this distribution it is In r rather than r which is normally 
distributed. An excellent discussion of this distribution is given by 
Kerker (Ref . 2) . 

The radius number distribution is given by 


n(r) 



exp 


/2tt p 3 r 


1 

2 


' lnr - In 
‘ p 3 



(8a) 
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The mode radius is 

2 

r m = p 2 ex P(-P 3 ) (8b) 

and its maximum value is 

P 1 2 

n(r ) ~ — exp(p„ /2) (8c) 

]/2v P 2 P 3 

The parameter P2 is the geometric mean of r and In p2 is the mean 
of In r. No series expansion could be found for the limiting. behavior 
of this distribution which tends rapidly to zero at both extremities. 

The log radius-number distribution is 


2,3 - Pi . 

_ i 

lnr - lnp 2 

2 ' 

exp 

~ o' 



/2tt p 3 

4 

p 3 

, . s 



It has a maximum at 


(8 d) 


r l B = P 2 ■ < 8e > 

so that 

V r lm ) = 2,3 Pl/(^ p 3 ) - (8 f) 

No series expansion could be found for the limiting behavior of this 
distribution which tends rapidly to zero at both extremities. 


The COSD is given by (Ref. 9, p. 183) 


N(r) ?= — erf c 


lnr - 1,: 


np. 


/2 P, 


(8g) 
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where erfc is the complementary error function. 

The limiting behavior of this function is as follows, 

As r -> 0, N(r) -*■ P x (8h) 

end, as r •> °°, N(r) -► 0 (8i) 

The parameter p 2 is the median for the COSD, i.e., 

N(p 2 )/N(0) = 0.5 

The moments for the distribution are given by (Ref. 8, No. 3*322/2) 

M k = P l P 2 k ex P (P 3 2 k 2 /2) (8 j ) 

The parameter p^ controls the polydispersity of the model and the 
parameter p^ has a multiplicative effect on mode radios. This model is 
graphically presented in Figs. (5A.1 - 5(3.3). 

Model 6 — The Normal Distribution (ND) Model 
The normal distribution is a symmetric distribution which is finite 
at r = 0 and, thus, strictly speaking cannot he used to represent aerosol 
SD's at small r. It can be used to represent SP's at other Ranges of r, and 
since it is a Gaussian distribution, which has well known properties, such 
a model aan be very useful in certain applications. It is given by 


n(r ) = 


/ 2tt 


r 1 
exp t - J 


r- 


£2. 


U P 3 


(9a) 
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Its mode radius is given by 


r m ^ P 2 


(9b) 


and the maximum value is 


n(r ) = — 

/2? P, 


(9c) 


Pi i p 0 

As r + 0 , n(r) a — exp( - i (-=■) 

v^rTp, 2 P 3 


2 ) 


j. (1 + rp 2 /p 3 ) (9d) 


No series expansion could be found for the asymptotic behavior of this 
function which tends rapidly to zero. 

The log radius-number distribution is given by 


n (r) 

JU 


2.3 r 

/zT p. 


exp 


' X. 
2 


r r -p 2 i 


3 j 


(9e) 


Its mode radius is given by 


r im = 


P 2 + ^P 2 2 + 4 p 3 2 


2 


( 9 f ) 


and the maximum value is 


2.3 p 


n (r 1 ) = — . 

L lm /Jir P, 


1 


nr. 2 

cP 2 + /p 2 +4p 3 


exp 


' 

2 , 2 ” 


'-1/2 

P 2 +4p 3 ' 


2 



L P 3 J 



(9g) 
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The limiting behavior of the distribution is as follows 

2.3 p^ r 

As r -0,n L <r)* jWJT 


(9h) 


and, as 'r -*■<», n T (r) -»• 0 

L 


The COSP is given by (Ref. 9 , p. 183) 


(9i) 


N(r) = — erfc 


r-p. 


/2 P J 


<9j) 


The limiting behavior of the distribution is as follows (Ref. 9, 
p. 183), 

P, P)' r-p 

As r -v 0 , N(r) .2 — [l -/ - ( - -■—)] (9k) 

2 ■ IT P3 


P 1 P 3 


and, as r 00 , N(r) ~ yjf' ( r _ p ) ex P 


r r-p. 


/^P 3 J J 


( 9 1 ) 


The moments of the distribution are given by (Ref. 8, No. 3.462/1), 


M k = m k! 6XP <_p 2 2/4p 3 2 ) D -k+l (P 2^ 21 


(9m) 


where D j. is a parabolic cylinder function. 

In the normal distribution, the parameter p 2 controls the mode radius 
and the parameter p^ controls the polydispersity . The model is graphically 
presented in Figs. (6A.1 - fjC.3). 



22 


Model 7 — The Generalized Distribution Function (GDF) Model 
This distribution is finite at r = 0 and thus does not, strictly 
speaking, represent particle size distrdbuitons at small r.„ However, it is 
a versatile function with a wide variety of applications, including altitude 
distributions, (Refs. 3 and b) and is therefore included here as a potential 
representation of aerosol SD's. 

The radius-number distribution is given by 

(10a) 


(10b) 
(10c) 

The limiting behavior of the distribution is as follows 

As r 0, n(r) - p (l+ r (p -l)/ p ( p +1)^ (lOd) 

1 2. 3 4 ■ J v' 

o 

and,. as ■ r ■* «> , n(r) - 0 as p.O+p.r exp(- r/p, ) " "(lOe) 

The log radius-number distribution is given by 

2 2 
n L d) = 2.3 p 1 (l+p 2 ) r exp(r/p 3 ) / (.p 2 * exp(r/p )) 


n(r) 


pj ( 1 + P 2 ) ex P (r/p 3 ) 


{ p 2 + exp (r/p 3 ) } 


Its mode radius is given by 


\ = p 3 ln »2 


and . the maximum is 


P, (1 + PJ 

n ( r ) = -i — £_ 

m / 

4 P 2 




(lOf ) 


Its, mode radius r lm is given by the solution of the equation 


p 2 a + r l» / « , 3 > * < r ln, /p 3 ' l) exp (r lm /P 3 ) 

apd the maximum is given by 


n L <r ln ,) = 2.3 Pl (Up 2 ) 2 <l lm -P 3 ) <l 1 „+P 3 ) / « P 2 r lm ) 

The limiting behavior of the distribution is as follows 
As r ->0, n(r) * 2.3 p x ( 1+P 2 ) 2 r 


apd> as r oo, n (r) ^ 0 as 2.3 p x (l+p 2 > 2 exp(-r/p 3 ) 


The COSD is given by 


Nit) « 


P 1 ( 1 -P 2 > 2 P 3 

1 1 1 I T — ; T-T" 

1 * 

P 2 + exp(r/p 3 ) 


The limiting behavior of the distribution is as follows 

As r -»• 0, N(r) * P 1 (l+p 2 ) p 3 (l-r/p 3 (p 2 +l)) 

and, as r + », N(r) -> 0 as p^ (,1+p ) ? exp(-r/p 3 ) 

The analytic expression fop the moments of the distribution 
cpuld pot be evaluated. 
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: (log) 

(10h) 

(lOi) 

(10j) 

(10k) 

dost,) 

( 10m) 
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The parameter can be considered as a scale radius and the parameter 
P 2 determines the type of function. For p 2 = 0 the distribution becomes 
an exponential, and for small p 2 the function initially falls off more 
slowly than the exponential. As the parameter p^ increases, the spread 
or polydispersity of the function increases. The model is graphically 
presented in Figs (7A.1 - 7C.3). 


Model 8 — Power Law Generalized Distribution' Function (PLGDF) Model 
This modal is a versatile function which is most useful when the data 
to be fitted have broad peaks. The radius-number distribution is given by 


n(r) 


P L exp 


(P 2 / r 


P4 +1 P 4 

r (l+p 3 (exp (p 2 /r ) 



(11a) 


Its mode radius is given by the solution of 


(1-P 3 ) (p 2 P 4 + (P 4 +1 ) r m ) = p 3 . exp (p 2 /r m ) (p 2 p 4 -(p 4+ l)r m ) 


(lib) 


and the maximum is 


n(r ) = 
m 


Pl < P 2 P 4 * (p 4 +1) r m > (p 2 p 4 ' (p 4 +1> r » > 

2 2 P4 + ^ 

4 P 2 P 3 (l-P 3 ) P 4 V 


(11c) 


The limiting behavior of the distribution is as follows 


P, 


As r -* 0, n(r) -v 0 as p^ exp (- p 2 /r ) 


(lid) 


and, as r -> e°, n(r) ~ p^ 


-(p 4 +D 


(lie) 


25 


The log radius-number SD is given by 


n L (r) = n 

L p 


2.3 exp (p 2 /r ) 


.4 


X (1+P 3 (exp(p 2 /r ) -1))‘ 


Its mode radius is given by the solution of 


Pv 


(1 -P 3 ) ( P2 +r lm } = P 3 ex P ( P 2 /r lm > ( P 2 ^ r lm } 


and its maximum is given by 


% 


W = 


P l (p 2 + 'l» > > 


4 p 2 P 3 U-P 3 ) r lis 

The limiting behavior of the distribution is as follows 

p. 

As ? G , n ]L (r) * 0 as 2.3 ?1 exp C-p 2 /r 4 ) 


and , as x w, n L (r) ~ 2.3 p^^ r 


The COSD is given by 


"P> 


N(r) 


H 

Pj^ (exp (p 2 /r ) - 1) 

l - ' 1 1 ' i "V -rr-" . . .r r--,,, - 

P 4 

P 2 P 4 (1+p 3 (exp(p 2 /r ) -1)) 


The limiting behavior of the distribution is as follows 


Po/r" 


As r' 0, N(r)'* P 1 /p 2 P 3 P 4 (l+l/p 3 (e -1)) 


and, as r -*■ oo ? N(r) 0 as — — ’ J 


P 2 p 4 


(exp(p 2 /r )-l) 


( l+p 3 <exp(p 2 /r ) —1 ) J 


(Ilf) 

dig) 

(llh) 

(Hi) 

dip 

(Ilk) 

du) 

(11m) 
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The moments, other than the zeroeth moment (N(0)) , cannot be 
calculated for this distribution. However, there is a special form of 
the function, referred to as Model 8B, for which higher moments can be 
calculated, but it has no analytic form for the COSD. The radius-number 
§D for this function is given by 


n(r) .= 


P L exp (p 2 /r ) 


r 4 ( l+pg (exp(p 2 /r 2 )-!)) 2 


(12a) 


Its mode radius is obtained from the, solution of the equation 


(l-p 3 ) (2p 2+ p 4 r m 2 ) =' p 3 C2p 2 -p 4 r m 2 ) exp(p 2 /r m 2 ) 


(12b) 


and its maximum is given by 


n(r ) = 

m 


P 1 (2 P 2 + P 4 r2)(2 P 2 ~ P 4 r m 2) 

2/i \ P4 

16p 2 p 3 (1 - p 3 )r m ... 


(l 2 c) 


The limiting behavior of this distribution is as follows 


As 


r 0, n(r) ->-0 as p^^ exp(-p 2 /r 4 ) 


and, as r °°, n(r) ~ r 


-P/ 


( 1 2d ) 

( 12e) 


The log radius-number distribution is, giyen, by 


n T (r) = ■ 

L p ,-l 


2.3 p L exp (p 2 / r ) 


r 4 ^ 1+ P 3 (exp(p 2 /r 2 ) - l)) 2 


( 1 2 f ) 


Its mode radius is given by the solution of the equation 
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(l-P 3 )( 2 P2 + (p 4 - 1 ) r ^ i ) = P3 (2 P2 _(p 4^ 1)r l m ) ex P(P 2 /r iV ) 


( 1 2g ) 


and its maximum is given by 




2,3 [, 1 (2p 2 t(p 4 -l)r 1 2 a )(2p 2 -(p 4 -l)r 2 lm ) 


16 


2 


p 2 P 3 a-P 3 )r 


P 4- 1 


lm 


The limiting behavior of the distribution is as follows 


As r -*■ 


2 

0 , n^(r) +■ 0 as 2.3 exp(.-p 2 /r ) 


and 


> as r -* oo f n L (r) ~ 2.3 pj^ p 


-(P 4 -1) 


( 12 j ) 


The moipents for this distribution are given by (Ref. 13, No. 313.11) 


M k ~ T~2 P4k 


2p 3 p 2 


r(p 4k ) ^ 

i =0 


M 


1 ~ 


p 3 


" (p 4k‘' 1) 

(i + 1 ) ( 12 k) 


p^-(k+l) 

where p 4k = — ^ , P 4 > k + 1 


The parameter behavior for both functions (Eqs. (11a) and (12a)) is 
similar. The parameter p 2 controls the rate of fall-off at small radii 
whilp the parameter p 4 controls the rate of fall-off at large radii. The 
parameter p^ controls the spread of the distribution, the breadth of the 
peaks increasing with large values of p^» The model is graphically pre- 
sented in Figs. (8A.1 - 8 F, 3). 
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GRAPHICAL CATALOG OF SIZE DISTRIBUTION MODELS 
Description of Model Catalog 

The model catalog consists of log-log graphs (6 x 4 cycles), and/or 
in some cases, semi- log graphs, on which n(r), n^Cr) and N(r) are separately 
plotted against radius (ym). For each of the aforementioned models, sets 
of plots have been produced, each set depicting the behavior of the model 
as a function of only one of its parameters, the other parameters being 
kept constant. Thus, each catalog graph (CG) contains several plots, one 
for each value of the parameter being varied, the parameter values being 
chosen to cover the range of variation in aerosol size distributions likely 
to be encountered. A CG with a p^, Pg or p^ variation will be hereafter 
referred to as a Pg CG, pg CG or p^ CG, respectively. Each catalog plot 
(CP) is labeled by a different symbol against which are printed the value(s) 
of the parameter (s) , outside the right-hand boundary of the graph. Also 
printed there, are the model number and the mode radius (RM) . For the 
first plot all the parameters are listed; and for subsequent plots only 
the values of the varying parameter and RM are printed. In cases where the 
SD has no mode radius, the value RM has been set equal to 0.01. The 
annotation on logarithmic axes gives the power of 10. 

For each set of parameter values, three corresponding CG*s for n(jr), 
n (r) and N(r) are presented. The values of the parameters held constant 
in each case, in general, were selected so that they were near the middle 
of their likely range. In all cases, the parameter p^ is a scaling factor 
which has been chosen so that the maximum values of n(r) , n^(.r), or N(r) 
are unity. This allows all of them to fit on one graph and facilitates 
estimation of parameters. 
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Use Of Model Catalog 

In this section* :a step-by-step procedure will be described illustrating 
how the catalog can be used to obtain the best-fit values of the model 
parameters. 

Step 1 : In case the two transparencies of log-log and semi-log graphs 

provided in the pocket on the back cover of this report are missing, take 

the blank log-log or semi-log graph provided at the end of the catalog and 
get some sharp photographic transparencies made of it. Xerox or vu-graph 
transparencies will not do as. they generally distort the graph unevenly. 

Step 2 : Next, plot the experimental data for n(r) , iw(r) or N(r) on 

the transparency, which. is ref eired to as the data transparency (DT) . If 
the data is not too noisy, one can visually detect some geometrical trend 
in the data points. A free-hand- smooth line may be drawn through the data 
points to accentuate the shape of the geometrical curve or line. 

Step 3 : The DT should then be overlayed on the different CG's belonging 

to the selected analytical model until a visual matching of the data points 

to one of the plots is obtained. With some experience, such a match can 

quickly be obtained. However, one often finds that the visual fit of the 
data points to a CP is reasonably close but not exact. The values of the 
parameters corresponding to the closest plot are then read off the CG, and 
used as initial parameter estimates in a nonlinear .least squares (NLLS) 
computer program or some other optimization code to obtain the best-fit 
values of the parameters. 

Since all the CP's for n(r), n_ (r) and N(r) are normalized to unity 

L 

at their individual maxima, the experimental data points may not, in general, 
fall in the same range as that of the CP's; thus vertical translation of the 
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DT will be necessary to obtain the visual match. In addition , horizontal 
translations of the DT, where permissible, may sometimes be necessary. In 
any case, care should be taken to ensure that no rotation pccurs between 
the CQ end the DT, 

As a rule, t^e horizontal translations are invalid when the parameter 
being varied controls the mode radius r^ and are valid if it controls the 
pplydispersity op the rate of fall-off at the extremities, 

PARAMETER ESTIMATION 

Therp are tlfree methods by which one can obtain, with relative ease, 
estimates of the parameters of the size distribution models that are 
close to the bestrfit values. In each case, the experimental data is 
plotted on a log- log pr semi-log graph paper. The first is based solely 
on the interpretation of the plotted data in terms of the mathematical 
properties (such as limiting behavior, mode radius, moments, etc.) of 
an analytic mpdel; the second is based solely on the visual matching of 
the plotted data with a set of parameterized plots of that model; and the 
third, which is the most versatile, is simply a composite of the first two 
methpds. It is the last methods, which will be discussed in this section. 
In this method, it is essential that the experimental data he plotted on 
the photographic transparency of the unused log- log graph provided with 
the catalog. In some cases it piay be helpfyl to use the semi-log graph, 
also provided with the catalog, 

A straight line on a log-log graph implies a power law, and on a 
semiylog fpaph, an exponential function, If for large r, the data points 
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on a log-log graph fall on a straight line, then the size distribution 

_P 2 

models which should be considered are those which vary as, say r for 

r -v oo, the slope being the value of the exponent, p ^ . This would immediately 

suggest that models 1, 2 , and 4 are the appropriate ones to work with. If, 

on the other hand, for large r, the data points on a semi-log graph fall on 

a straight line, the models which should be considered are those which vary 
-P 2 r 

as, say e as r-»- ®> , the slope being the value of the coefficient 
This would suggest Model 3 with p^ =1 or Model 7. If, however, the 
experimental data do not fall on a straight line when plotted on either 
the log-log or the semi-log graph, then Models 3 and 5 are important. In 
the following discussion on parameter estimation, primed symbols will 
refer to the data estimates and unprimed symbols, to the CP values. 

The parameter pj is always the last parameter to be determined. 

Since the estimation of the scaling parameter p| is essentially the same 
for all models, it will be discussed before describing the procedures for 
estimating other parameters for each of the models. The estimates for pj 
can be obtained by one of the following three methods, after the estimates 
for p^, p^ , etc. have been obtained. 

(a) If one assumes a linear relationship between the observed ( y °. ) 
and calculated (y^) values, i.e., 

y° = p{ y': 0.3a) 

J L J 

then by using linear regression methods it can be shown that 


P] =. ( X 

j = l 



(13b) 
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where J = number of data points, y ^ = observed SD value at r 
c 

7 ^ ? calculated §D value at r ^ , using pj = 1.0 and the p^, p^ 
already determined. 


(b) 


(observed value at r . /calculated value at r.) 


and 
etc . 


(13c) 


where the calculated value of the SD is obtained by substituting 

pj - 1.0 and p^ f ? etc. as estimated from the CP's, in the expression 

for the SD. If possible, the values at the maximum should be used. 


(c) pj = (observed value at r^/catalog value at r^.) x p^ 


(13d) 


Evep though the second method might perhaps be the easiest for 
estimating p J , the first method is found to be more convenient when a 
computer optimization code is used, since the formula in Eq. (13b) can 
easily be incorporated into the program. The third method seems to be 
tf>e most complicated since often several CP's with different p^ values 
are used to estimate the parameters, P 2 / P 3 r etc. 

The first step is to determine the other parameters p^j P 3 , etc. 
for each of the models. The procedures for estimating these parameters 
are described under the subheadings for the different models as follows. 

Model Is The slope of the straight line on a n(r), n. (r) or N(r) 

' _r "' ' ' ' Ju 

plot determines the parameter P^* 

Model 2 : This model behaves as a power law for both large and small 
r. Therefore, in general, it is best to determine p^ and p^ first, and 
then determine p,J. The parameter estimation differs somewhat for the n(r), 
n^(r) or N(r) data plots. 

(a) n(r) data, pj can be obtained from the positive slope 
(= P 3 ~ 1) of the data plot; p^, from the negative slope (= -1 - p^ (p^ -1)), 
after substituting for p^ ; and p^ is found by substituting for r^, p^ 
apd p^ in EJq, (5b). 
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To obtain pj using the catalog, the DT is overlayed on 
the p^- CG (Fig. 2A. 3) , and translated along x- and/or y- axes until 
the data points corresponding to the small limit of r match, as closely 
as possible , one of the CP's. 

To obtain p^, the DT is overlaid on the p^ CG (Fig. 2A.4) 
and translated parallel, to th,e x- and/or y- axes until the data points 
for large r closely match one of the plots. (Since both p^ and p^ ; 
affect the slope for r , it may be necessary to adjust the p^ value,, 
using, . 

P 3 (P 4 “!) = P 3 (P 4 -U (14a) 

if the p^ value for the CP differs from p^.) Having obtained p^ and 
p^ , estimates of p£ can be obtained from Eq. (5b) or from the CP's. In 
the latter case, overlay the DT on the CG for which (p^ , p^) are as 
close as possible to (p^, P 4 ) (Fig. 2A.1 or 2A.2) so that the data 
maximum coincides with 1.0 on the CG y- scale (CP's are normalized 
so that their maximum value is 1.0). Translation of DT along the x-axis 
is not permitted, since the parameter being estimated controls r m . If the 
data peak does not coincide with any of the CP peaks, but lies between two 
CP peaks then p^ can be obtained by the following interpolation relations. 


p'/r ' = p./r 

2 m 2 m 


(14b) 


(If p^ and/or p^ for the CG differ from those values for the DT, then use 

1/Po 1/Po 


p 3 " 1 1 

— t 

p 3 - 1 | 

1 + p^(p^-l) ^ 

- P 2 

1 + p 3 (p 4 -1) 


( 14c) 


to obtain a better estimate, p^', for P 2 >) 
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(b) n^ir) data. The n L (r) data is plotted on the log-log 
transparency. The determination of the parameters is done in essen- 
tially the same manner as for the n(r) data. As r 0 , = positive 

slope; as r -»■ 00 , = 1 + negative slope/p^J and p£ is found by sub- 
stituting for rl , p^ and in Eq, (5g) . The parameters p^ and p^ 

are determined from Figs- (2B.3 and 2B.4) respectively, in the same 
manner as for n(r). Similarly, p^ is obtained in the manner described 

for n(r) data using Figs. (2B.1 or 2B.2), provided r and r' are re- 

mm 

placed by r^ m and r^ in Eq. (14b) and (14c) is replaced by 


-i/p; - 1 /P 3 


P 4- 11 

> p 2 

f ^ 

P 4- 1 

J 


. 


(14d) 


(c) N(r) data. The parameter determination from the N(r) 
dqta is more difficult than from the n(r) data. Plot the N(r) data on the 
log-log transparency, visually fit a line through data points, and, 
if possible, extrapolate it to r = .Ollim. 

To obtain p^ , overlay the DT on the C 0 (Fig. 2C.1 

or 2C.2) whose r -»■ <» behavior best matches that of the data, so that 
the N(0) for the data coincides with the N(0) for the CP. If data 
lies between fwo CP'Sj then qse the following extrapolation formula 
to obtain p£, viz. 



(14e) 


where r and r’ are radii with the same normalized value N(r) on the 
straight portion of the curve (Fig. 9) and s' and s are the slopes of 
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the straight portion of the curve. It is best to use the CP which turns 

closest to the data. P! can then be obtained from the relation 

4 



( 14f ) 


To estimate overlay the DT on the p^~ CG (Fig. 2C.3) and translate 
along x- and y-axes to obtain the closest match between data points 
for r 00 and one of the CP's. (If p^ for the CG differs from p^, then 
to obtain a better estimate for p^, use Eq. (14a).') 


Model 3; In representing the size distribution data with Model 3, 
one generally obtains estimates of the parameters p£ and p^ first and 
those of p^ from them. 

(a) n(r) data: If n(r) data is available for r 0 , then 

one can obtain p^ from the positive slope (=p£). To obtain p^, overlay 
the DT on the p^~ CG (Fig. 3A.1), and translate it along x- and y- axes 
until a CP matches with the slope of the data for r ->-0, where the positive 
slope = p . To obtain p ' , overlay the DT on the p -CG (Fig. 3A.4) , and 
translate it along x- and y- axes , until a CP matches with the data for 
r^- 00 , for which limiting behavior is an exponential type fall-off. 

P^ can then be obtained by substituting r^, p^ and pi, in Eq. (6b) or from 
the appropriate CP. In the latter case, overlay the DT on a CG with (p£> p^) 
closest to (p^, p^) (Fig. 3A.2, or 3A.3), so that the data maximum 
coincides with 1.0 on the CG y- scale. The following interpolation formula 
can be used to obtain the p^ if the data do not coincide with one of the 


CP peaks. 
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, P 4 


m 


m 


(15a) 


(If p 2 ^ and/or p^ i p^, then 


, , 1/p 4 


-’3 >; 


i/p, 


p 3 P 4 


(15b) 


Qian be used to pbpain a better estimate p!J for the p 3 parameter since p£ and 

p! also influence the mode radius.) 
r 4 

(b) n^Cf) data: Determination of the iiodel 3 parameters 

from, the n T (r) data is similar to that described for the n(r) data, 
provided (p£ + I) is used for the positive slope, p^ is estimated from 
p£ and p^ using the expression for r^ given in Eq. (6g) , r m and r^ are 
replaced by r^ m and r£ m in Eq. (15a), and Eq. (15b) is replaced by the 
relation 


p 2 + 1 

*T*T 


1/p 4 


1/p, 


P 2 + 1 
P 3 P 4 


(15c) 


(c) N(r) data: Determination of the model parameters frpm 

the N(r) datp is difficult. It will generally be necessary to extra^ 
polate the data to r = Q.Olym. 

To obtain p 2 overlay the DT on a CG with plots whose 
Slopes for r -> °°(Fig, 3C.1 or 3C.2) are similar to that for the data, 
so t;hat^ the extrapolated N(0) value for the data coincides with N(0) 
for the CP's. Then p^ is the p 2 ~ value for the CP which best fits the 
data particularly with regard to the curvature. It may be necessary 
to interpolate visually in order to obtain a better value for p 2 » 
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To obtain p^, overlay the DT on the p^- CC> (Fig. 3C.5) 
and translate it in a manner similar to that for the n(r) data, until 
a best fit to the data for r •* 00 is obtained. Similarly, p^ is obtained 
by overlaying the DT on the CG with (p£ , p^) closest to (p^j P^) (Fig* 
3C.3 or 3C.4), and translating vertically so that N(0) for the data and 
CPs coincide. Then determine which CP best fits the data, interpolation 
by eye may be necessary to give a better estimate of p^ • 

Model 4 : Determination of parameters for this model is similar 

to that for Model 3. In general, determine p£ and P 4 first, and then 

P 3 ’ 

(a). n(r) data: The parameter p£ can be determined from the 

slope of the data as r 00 . (slope = p£). To determine pj( from the 
CG's, overlay the DT on the CG (Fig., 4A.1) and translate 
the DT parallel to the x- and/or y— axes until the data points as 
r -coolie on or close to one of the CP's. The P 2 value for this CP 
then becomes the estimate p^« To determine p^, overlay the DT on the p^- 
CG (Fig. 4A.4) and translate the DT vertically and horizontally until 
the data points as r -* 0 lie on one of the CP's. 

The parameter p^ can then be determined by substituting 
for p^ , p^ and r^ in Eq. (7b) or using the catalog graph. To determine 
p^ from the catalog, overlay the DT on the p^- CG which best matches -the 
behavior of the data both as r + 0 and r -*■ °° (Fig. 4A.2 or 4A.3) so that 
the data peak coincides with 1.0 on the CG y- scale. If interpolation 
between CP values is necessary, use the following relationship: 
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P 3 


r' 

m 

r 

m • 


, P 4 


(16a) 


(If the values of p 2 and/or for the CG differ from the estimates p£ and 
p^, then a better estimate p'^ can be obtained by using the relationship: 


(16b) 


- 

i/ pi 

1 

' 

P 3 P 4 


P 3 P 4 

p;' 
v 2 


l P 2 


where p 2 is the new estimate for p^.) 

(b) n (r) data: This is very similar to the parameter deter- 

v ■ 

minatiop fqr n(r). The slope as r -* “can be found from the data (slope 
s P 2 rl). Using the catalog, p 2 and are determined from Figs. 4B.1 
and 4B.4 respectively in the same way as p 2 and p^ for n(r). The 
parameter p^ is also determined in the same manner as for n(r) however, 
if interpolation between CP's is needed, r^ m and r| m should be used 
instead of r^ and r^ in Eq. (16a). (If (p 2 , p^) for the CG differ from 
the estimates (p 2 , p^), a better p^ estimate can be obtained by using 

1/P/, 


i/p i 




f P 3 P 4l 

k- x J 


P 2 _1 


(16c) 


where p” is defined as the npw estimate for p^.) 

(c) N(r) data: Determination of parameters is generally more 

difficult than for n(r) or n^(r). The slope as r 00 is given by - 
(p^-1). General ly , it will be necessary to extrapolate the data to 
r = O.Olym. 
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To determine from the catalog, overlay the DT on the 
P 2 ~ CG (Fig. 4C.1) and translate the DT vertically and horizontally until, 
the best fit CP for the region r -> 00 is found. To determine p^ , overlay 
the DT on the p^- CG (Fig. 4C.4) so that N(0) for the data coincides with 
N(0) for the CP's. The estimate for p^ will be the p^ value for the best 
CP, giving particular attention to the fit in the region of curvature. 

To determine p^ , overlay the DT on the p^- CG which best approximates the 
behavior of the data as r -> «> (Fig. 4C.2 or 4C.3) so that N(0) for the 
data coincides with N(0) for the CP's. The estimate for p^ is the value of 
Pg for the CP* which best fits the data. The point where the curve turns 
is of most importance here and some interpolation by eye may be necessary. 

Model 5 ; In representing the size distribution data with Model 5, 
it is best to obtain an estimate of p^ first and use it to determine p£. 

(a) n(r) data: If sufficient data are available , both p,' and 

p^ can be determined from the mathematical properties as follows : 

P 3 = 0.424 ln(r 1 /r 2 ), ^ > r 2 (17a) 

where r^ and r 2 are the radii at ha If -maximum points. Then p^ can be 
obtained by substituting for p^ and in Eq. ( 8 b). To determine the 
parameter using' the catalog, overlay the DT on the p^- CG (Fig. 5A.3) 
and translate it so that the maxima for the data and a CP coincide. 

Then p^ is the same as the p^ value for the plot which best fits the 
data. Attention should be paid to both the spread of the data, and 
the rate of fall-off as r °°. Often some interpolation may be 


necessary 
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Once pj is obtained, can be determined by substitution 
in Eq • (3b) or from the catalog. In the latter case, overlay the DT on 
the Pg*- CG wifih p^ closest to p^ (Fig. 5A.1 or 5A.2) so that the data peak 
coincides with 1.0 on the y- scale of the CG. If the data do not coincide 
with one of the CP peaks > p£ can be estimated from the closest one by using 

r' 
m 

r (17b) 

m 

(If p^ 4 p^> then a better estimate for p 2 can be obtained from 

p^ e,xp(-pj 2 ) = p 2 exp(-p 2 ) (17c) 

where pi,' is the new estimate of P 2 for the data.) 

(b) n ^(r) data: can easily be obtained from the data peak, 

since p 2 = r^. Determination of the model parameters from the CP's is 

essentially the same as for the n(r) data, provided r and r' are 

m m 

replaced by r^ m and rj m in Eq. (17b). Note that no adjustment of p^ 

Kill be necessary if p^ for the catalog differs from p^. 

(c) N(r) data: Generally, it is best to extrapolate the N(r) 

data to r = O.Qlytn to obtain an estimate of N(0). Then p£ can be 
determined from the fact that p 2 is the median, i.e., N(p 2 ) = 0.5, or 
in the case of unnormalized data, N(p 2 )/N(0) = 0.5. 

To obtain p^ , overlay the DT on the p^- CG (Fig. 5C*3), 
translate it horizontally and vertically until a match of the data for 
r -> oowjth a CP is obtained. (If p^ has been estimated as described 
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earlier, then the data point ( p ^ , 0.5 N(0)) should coincide with the ! 
point ( p 2 > 0.5) on the CP. If the median points for the data and the _ 
CP do not coincide it probably means that the original estimate for p£ 
is not very good and a better one would be the radius on the DT which 
corresponds to the median point for the CP.) 

To determine p^ using the catalog, overlay the DT on the 
P 2 ~ CG with Pg closest to p^ (Fig* 5C.1 or 5C.2) so that the N(0) values 
for the data and CP ' s coincide . is then the P 2 ~ value of the . CP with 

which data points match. Interpolation by eye may be necessary. 

Model 6 : Since this model was included in the catalog for complete 

ness only, a description of how to determine parameters from the catalog 
will not be given. However, if the model is to be used, the data should 
be plotted on a semi-log graph and the procedure for determining the 
parameters will be similar to those described earlier. 

Model 7 : Data for this model should also be plotted on a semi-log 
graph. As in the case of Model 6, this model was included for the sake 
of completeness and will not be discussed here. However, determination 
of the parameters for this model is discussed by Green ; (Ref. 6). 

Model 8 : Two forms of this model were discussed earlier, but para- 

meter determination is the same for both, so that separate discussions 
will not be necessary here. 

(a) n(r) data: To determine p£ , overlay the DT on the P 2 ~ 

CG (Fig. 8A.1 or 8D.1) so that the radii are aligned and the data peak 
coincides with 1.0 on the CG y— scale. p^ is the same as P 2 for the 
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CP which coincides with the data peak. Interpolation by eye may be 
necessary to get a better estimate of p^. 

To determine p^ the DT is overlayed on the p^~ CG 
(Fig. 8A.2 or 8D,2) and translated vertically and horizontally until 
the CP which best matches the peak width is found. It is this para- 
meter which controls the broadness of the peak and is of importance for 
data with broad flat peaks. The parameter p^ can be found from the slope 
of the data for r ->■ 00 ( = -(p^ +D for Model 8, = -p^ for Model 8B) or 
from the catalog. Using the catalog, overlay the DT on the p^- CG 
(Fig. 8A*3 or 8D.3) and translate vertically and horizontally until a 
CP which matches the data for r -> <» is found. 

(b) n^(r) data: The procedure for determining p^, pg, * s 

exactly the same as for n(r) except that in the case of Model 8 Figs. 8B,'l 
8B,3 are used and the slope for r -> °° is -p^ and in the case of Model 8B 
Figs. 8E.1 - 8E.3 are used and the slope for r -* 00 is and -p^ - 1. 

(c) H(r) data: Extrapolate the data to r = O.Ojym, then 

to determine p£ , overlay the DT on the p£- CG so that the radii are 
aligned and N(0) for the data coincides with N(0) for the CP's. The 
estimate for p^ is given by the value of p ^ for the CP which best 
matches the data, particularly in the region where it turns. 

To determine p^ overlay the DT on the p^~ CG (Fig. 8C.2 
or 8F.2) so that N(0) for the data and the CP ' s coincide . The estimate 
of p^ is the value of Pg for the CP which best matches the spread or 
polydispersity of the data. To determine p^, overlay the DT on the p^- 
CG (Fig. 8C.3 or 8F.3) and translate the DT vertically and horizontally 
until the GP which best matches the data for r '->■ 00 is found. 
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Multimodal Size Distribution Models 

Quite often the SD data may indicate that the distribution is multi- 
modal (Refs. 14 and 15). Such SD's present no particular difficulty in 
representing them analytically in terms of the aforementioned SD 
models. One can represent multimodal size distributions simply by 
adding appropriate mathematical functions, each component term repre- 
senting a peak in the size distribution data, such as, 

J 

n(r) = l a. n (r) (18a) 

j=l J J 

where J is the number of component terms and the a^ represent the 
adjustable constants and are equivalent to p^ for the unimodal models. 

In describing the procedure for fitting multimodal data using the 
catalog, a bimodal model will be condensed for simplicity. 

Step 1: Plot the data on a log-log transparency and dot in the 

probable behavior of the two modes as illustrated in Fig. (10). 

Step 2: Decide the most appropriate model for each mode using 

the guidelines already outlined for the unimodal models. 

Step 3 : Determine the parameters for each mode using the guide- 

lines outlined earlier for unimodal models. 

Step 4: Determine the scaling factors for each mode. The simplest 

method for determining the scaling factors is the solution of the 
simultaneous equations 

n ( r 1 ) = a^^rp + a 2 n 2 (r 1 ) (18b) 

n(r 2 } = a l n l (r 2 ) + a 2 n 2 (r 2 } 


(18c) 
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where a 2 are the scaling factors required, and ^ and r^ are the 
mode radii for modes 1 and 2, respectively.. 

In some cases, it may be difficult to get good estimates for a^ and 
a 2 j either because of uncertainty in the accuracy of some estimates or 
because the scaling factors differ by several orders of magnitude. An 
alternative approach to determining the scaling factors under these 
circumstances would be to assume that each mode makes a negligible con- 
tribution to n(r) in the region where the other mode dominates and, thus, 
the qcaling factors can be found by using the following relationships: 

n(r^) = a^njXr^) (18d) 

n( r 2 ) = a 2 n 2 (r 2 ) (18 e ) 

whqre r^ and r 2 are radii at which modes 1 and 2, respectively, dominate. 

Examples of Analytic Representation of SD Data 
Tropospheric Aerosol Size Distributions (Junge Data) : Junge ' s 

data (Ref. 1) for the average size distribution n(r) of continental 
aerosols in the altitude region (0-3 km) were plotted on the log- 
log graph transparency. One could readily see that for large radius 
particles, the n(r) data points lay along a straight line and suggest 
thereby a power law behavior. Thus, Models 2, 4, and 8 seemed 
appropriate. In the following, examples are given for fitting the SD 
dafa with these three models. Model 5 was also fitted to the data to 
illustrate the poor fit which can be obtained if an inappropriate 
model is chosen. 

Model 2: According to our earlier prescription, the DT was 

overlaid on the p^ - CG (Fig. 2A.4), and translated vertically 
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_q _2 

until the point (10 , 1.0) on the DT coincided with (10 , 0.1) on 

the CG. The data points for large radii fell on the p^ = 2.0 plot; 

thus, Pg = 3.0 and p^ = 2.0 (the CP parameters) represent the data 

for large radii. Since there is not sufficient data at the smaller 

end to determine p^ from the slope in that region, concentration will 

be on the curvature in the region of the peak. 

Still using Fig. 2A.4, the DT was translated vertically 

6 

and horizontally so that the point (0.5, 10 ) on the DT coincided with 
(100, 1.0) on the CG. Although p^ = 3.0, p^ = 2.0 represented the large 
r data well, the CP had a sharper peak than the data points and thus 
suggested that p^ = 3.0 may be too high. 

The DT was then overlaid on the p^-CG' (Fig. 2A.3), and 

6 

translated vertically and horizontally so that the point (5.0, 10 ) on the 
DT coincided with the point (100, 1.0) on the CG. The plot for p^ = 2.0 
gave the best fit to the data. It did not fall off as sharply at small 
radii and the slope at large radii was a little too steep; however it 
did give a reasonable fit to the broad peak. 


The estimate of p^ had to be adjusted in the light of 
the new estimate for p^ = 2.0 since the slope at large radii depends 
on both p^ and p^. This was done by substitution in Eq. (14a) which 
gave p^ = 2.5. The only parameter left to be estimated was p^ which 

_3 

can be found from Eq . (5b). Since r T = 6.5 x 10 for these data, 

' . m 

substitution for r^, p^, and p^ in Eq. (5b) gives p^ = 0.013. The 
accuracy of this estimate was confirmed by estimating p^ from the 


catalog. The DT was overlaid on the P 2 ~CG (Fig. 2A.1), and trans- 

-2 5 

lated vertically so that the point (10 , 8 x 10 ) on the DT coincided 

-2 

with (10 , 1.0) on the CG. It was evident that the p£ for the data 
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was less than 0.05. The estimate for was found by substitution in 

Eq. (14b). Substituting p 2 (=0.05) and r m (=0.022) for closest catalog 

plot and (= 0.0065) for the data points, yielded p^ = 0.015. However, 

p 3 = 2'°’ p 4 = 3.0 on the CG and the estimates for the data were p* = 

3 

2.0, p^ = 2.5; thus, the p^ estimate obtained had to be adjusted by 

using Eq. (14c), which gave p 2 = 0.013, the same as the earlier estimate. 

These estimates (p* = 0.013, p^ = 2.0, p^ = 2i5) were 

used as the initial estimates for a nonlinear least squares (N LL S) 

fit. The results are given in Fig. (11) and the best fit estimates 

-2 

ate = 3.172 x 10 , p^ = 1.82, p^ = 2.65. It can be seen that 

the initial parameter estimates represent the SD slope and the curvature 

quite well. Note that the estimated curve has the right shape but is 

too far to the left of the data points. This is due to the broad peak, 

the estimate of pi was too low, since the actual r value was used. 

2 m 

It might have been more accurate to take r^ = mid point of peak region. 

The final best-fit values, obtained by the NLLS method, give a good 

analytic fit to the SD data. It should be noted that the greatest change 

occurred in p^, which effectively gave a higher r m value. 

Model 4 : Since the IMGD behaves as a power law for 

large 4, it is a suitable model for the Junge data (0-3 km). From the 

data points the slope as r + ® was found to be -4; thus the estimate for 

pi, was 4.0. To confirm this estimate, the DT was overlaid on the p 2 ~ 

-3 

CG (Fig. 4A.1) so that the DT point (10 , 1.0) coincided with the CG 

—2 —5 

point (2 x 10 ,10 ) , their respective x- and y-axes remaining mutually 
parallel. It was noted that the data points lay along the p 2 = 4.0 CP. 
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The DT was then overlaid on the p^- CG (Fig. 4A.4) so 

_3 

that the point (10 , 8.0) coincided with the CG point (0.3, 

_5 

10 ) . The data maximum then coincided with the maximum for the p^, = 

2.0 CP and the data points in the region of the peak fell on that 
curve. Thus, p^ = 2.0 was the estimate. 

Substituting the estimates for p^ and p^ and r^ — 6.5 

-3 -5 

x 10 in Eq. (7b) gave p^ =8.45 x 10 . To confirm the p^ estimate, 

the DT was overlayed on the p^- CG (Fig. 4A.3) so that the DT point 

(10 •, 8 x 10 + ^) coincided with the CG point (10 1.0). From this, 

it was evident p^ << 1.0. Substituting r^ for the data, and p^ and r 


m 


-5 


for the closest CP plot in Eq. (16a) gave p^ = 4.23 x 10 . However, 

(p 2 » P 4 ) f° r tbe CG were different from the (p^, P^) estimates for the 
data; thus, the final estimate for p^ was found by substitution in 
Eq. (16b), which gave p^ = 8.45 x 10 the same as the earlier estimate 
Figure 12 shows the results, the final estimates were p^ = 4.43, 


p^ = 0.788, p£ = 0.514. 

Model 5 : Since these data fall off as a power law for 

large r, the log-normal distribution is not the most appropriate model 
to use. However, it is used here in order to illustrate the poor 
fit one obtains when the model chosen is inappropriate and alsd to 
illustrate that such a model can sometimes be used to fit the data 
in a portion of the data range. 

The DT was overlaid on the p 3 ~ CG (Fig. 5A.3) and 
then translated horizontally and vertically so that the DT point 
(10 , 8.0) coincided with the CG point (6 x 10 10 "*) . At this 

point, the maxima for the data and the CP for p^ = 1.0 coincided. 

The data points lay almost exactly along the p^ =1.0 CP, except for 
very small r^ Thus, p^ =1.0 was the estimate. 
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The parameter was then estimated by using Eq. (8b) and 

_3 

r^ = 6.5 x 10 . Substitution in this equation gave p^ = 0.0176. To 

confirm this estimate from the catalog plots, the DT was overlaid on 

-2 5 

the p - CG (Fig. 5A.2) so that (10 , 8 x 10 ) on the DT coincided with 

_2 

(10 , 1.0) on the CG. It was noted that f° r t ^ ie data was less than 

0.1; thus, in order to obtain an estimate for an extrapolation was 
made by Eq. (17b). Substitution of the p^ and r^ values for the closest 
CP and r^ for the data gave p^, = 0.0176, the same as the previous estimate. 
Thus, the parameter estimated for this model were p^ - 0.0176, p^ = 1.0. 
Figure 13 shows the results of an NLLS fit using these 

-2 

parameter estimates. The best fit estimates were p^ = 1.208 x 10 
p^ = 1.21. The initial estimate provides a good fit to the data up to 
about 2ym, beyond which it falls off too steeply. The final best-fit 
values give a much better fit in the region 0.2ym - 30ym, but do not fit 
the peak region. As was mentioned earlier, the log normal model is not 
as suitable for these data as Models 2 and 4. However, if only a portion 
of the size range covered is of interest, say up to 2ym, or beyond 0.2 m, 
then this model is as good as any of the others. 

Model 8B : Since Model 8B behaves as a power law for 

large r it is a suitable model for the Junge data. It also has the 

advantage of being able to model broad peaks. To determine p^, the 

_2 

DT was overlaid on the p£- CG (Fig. 8D. 1) so that the point (10 , 

5 -2 

8 x 10 ) coincides with 10 , 1.0) on the CG. It was then noted 

-4 -5 

that the data points lay about halfway between the 10 and 10 CP's. 

The estimate for p^ was thus 5 x 10 

To determine p^, the DT was overlaid on the p^- CG 

-3 5 

(Fig. 8D.2) so that the point (5 x 10 , 8 x 10 ) coincided with the 
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point (2 x 10 , 1.0) on the CG. The data points coincided with the 

Pg = 40.0 CP; thus, the estimate was p^ =40.0. 

The parameter p^ was determined by overlaying the DT 

— 9 6 

on the P^- CG (Fig. 8D.3) so that the point (10 , 10 ) coincided 

-2 

with (2 x 10 , 10.0) on the CG. Since the data points lay along the 

p^ = 4.0 CP, the estimate for p^ was 4.0. 

Thus, the parameter estimates for this model were p^ = 

5 x 10 \ p^ = 40.0, p^ = 4.0. The results of the NLLS fit to this data, 

using these estimates, are given in Fig. 14. The best fit estimates 

were p^ = 4.700 x 10 p^ = 38.8, p^ = 4.03. As can be seen from Fig. 14 

the initial estimates represent the shape of the data quite well but the 

n(r) values are consistently too low. The final estimates give a very 

good fit to the data and model the broad peak very well. 

Cumulative Size Distribution : The SD data was obtained over 

Alaska at an altitude of 13 km, by Quartz Crystal Measurement 

(QCM) method . 1 The SD is in the form of the cumulative. 

Model 2 : The data points fell on a straight line 

as r -> «>; thus. Model 2 is an appropriate model. In this case, p^ 

was determined first, and then p^ and p^ which control the slope for 

large r. The first step was to extrapolate the data to r = 0.01. 

3 

This gave a value 8 x 10 for N(0) . 

The DT was overlaid ,on the p 2 ~ CG (Fig. 2C.1) so. that 
N(0) for the DT coincided with N(r) = 1.0 on the CG. It was noted 
that in this case the rate of fall-off of the data for large r was 

^ D . Woods, private communication. 
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the same as for the CP' s but that the data lay between the 3 x 10 
-2 

and 5 x 10 CP's. By using the interpolation formula in Eq. (14e) at 
-4 

N(r) = 3 x 10 (on. the CG) ,. assuming s = s' and substituting r = 

0.4, r' = 0.3 and p ^ - 0.05, the estimate p), = 0.0375 was obtained. 

To obtain an estimate for p^, the data values N(p^) and N(0) were 
substituted in Eq. (14f). This gave p^ = 4. 

The DT was then overlaid on the p^- CG (Fig. 2C.3) so 
that the point (0.5, 10^) coincided with the point (40, 10) on the 

CG. The data points lay on the p^ = 2.0 CP making p^ = 2.0 the estimate. 

However, both p^ and p^ affect the slope as r -*■ 00 , and p^ on the CG was 

different from the p^ estimated above, thus the estimate for p^ was 
adjusted by using Eq. (14a). Substitution gave p' ■ 1.33 as the final 
estimate. 

Thus, the parameter estimates were p^ = 0.0375, p^ = 1.33, 
p^ = 4.0. Figure 15 shows the results of an NLLS fit using the estimates 
for , p^ and p^ given above; the best estimates for the parameters were 
p), = 0.0434, p^ = 1.29 and p^ = 4.31. As can be seen from the plot, 
the initial estimate explains the data very well and the final estimates 
for the parameters are very similar to the original estimates. 

Model 3 : The DT was overlaid on the CG (Fig. 3C.1) 

_o 3 _2 

so that the point (10 , 8 x 10 ) coincided with the point (10 , 1.0) 

on the CG. The shape of the data curve in the region where it curves was 

very similar to the shape for the CP ' s , however , p£ < 1.0. In the absence 

of any reliable extrapolation formula, p^ = 1.0 was used as the best estimate. 
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The DT was then overlaid on the p^- CG (Fig. 3C.5) so 

that the point (10 , 10 -1 ) coincided with (3.0, 10 ) on the CG. The 

data lay close to the p^ = 0.3 CP so p^ = 0.3 was used. 

To obtain an estimate for p^, the DT was overlaid on 

-2 3 

the p^- CG (Fig. 3C.3) so that the point (10 , 8 x 10 ) coincided with 

_2 

the point (10 , 1.0) on the CG. The data points lay closest to the 

p^ = 30.0 CP; so this was taken as the estimate for p j . 

Thus, the estimates for the model are p^ = 1.0, p^ = 30.0, 

and p^ = 0.3. Figure 16 shows the results of an NLLS fit to the data 

using these estimates; the best estimates were p^ = 0 . 999 , 

P 3 = 30.0 j = 0.249. The initial estimate gave a reasonable fit to 

the data, and the final estimate fitted the data very well. 

Model 5 : When fitting cumulative data with the LND, it 

is generally best to extrapolate the data to obtain N(0), and use this 

3 

to calculate p£* Extrapolation of the data gave N(0) ~ 8 x 10 , and by 

3 

using N(p^)/N(0) = 0.5, it was found that N(r) = 4 x 10 at r = 0.02 and 

gave a first estimate for p^. 

The DT was then overlaid on the p^- CG (Fig. 5C.3) so 
-2 -3 

that (2 x 10 , 4 x 10 ) coincided with (1.0, 0.5) on the CG (i.e. the 

median points coincided). This suggested 0.5 < p^ < 1.0, since the data 
points lay between the CP'shaving those values. Looking at the behavior 
as r -*■ 00 , it was noted that the CP for p^ = 0.5 fell off much more 


sharply than the data. The DT was translated horizontally and vertically 


This Page Intentionally Left Blank 



53 


The DT was overlaid on the p^- CG (Fig. 3A.4) so that 

-2 2 -2 
the point (10 , 6 x 10 ) coincided with the point (3 x 10 , 10) on 

the CG. Then, by looking at the behavior as r + <», it was noted that the 

data points lay on a line which was approximately parallel to the p^ = 

1.0 CP. The DT was translated horizontally so that the point (20.0, 

3 2 

6 x 10 ) coincided with the point (10 , 10) on the CP. The data points 
did, in fact, lie close to the p^ = 1.0 curve and thus confirm the 
estimate p^ = 1.0. The parameter p^ was then estimated from r^, p^» 
and p^, by substituting into Eq. (6b); by taking r^ = 0.1275, this 
gave p^ = 78.4. 

To confirm the estimate for p^ the DT was overlaid on 

-2 2 

the Pg- CG (Fig. 3A. 2) so that the point (10 , 6 x 10 ) coincided with 

_2 

the point (10 , 10) on the CG. The maximum for the data occurred at 

the same r as that for the p^ = 10.0 CP; however, P 2 and p^ for the CG 

were different from the estimates p^ and p^ so that p^ estimate was 

adjusted by using Eq. (15b). Substituting P 2 , Pg> p^ for the CG and p^, 

p^ for the data gave p^ = 62.5. Thus, the original estimate for p^ 

seemed reasonable. Therefore, the estimates for this model were p^ = 

10.0,' Pg = 78.4, and p^ = 1.0. Figure 18 shows the best fit obtained; 

the final estimates were p^ = 8.38, p^ = 64.7, and p^ = 0.758. 

Bimodal Size Distributions : The data used in this example 

were collected as part of an investigation of high temperature 

combustion aerosols (Ref. 17). The original data were in the form 

of the aerosol mass in each of eight diameter groups, ranging from 

0.43 ym to 20 ym. For this example, the masses have been converted 

to number of particles by assuming spherical particles with density 
3 

1 yg/cm and average diameter equal to the midpoint of the diameter 
group in question. 
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Figure 19 shows the original data, with a probably curve for 
the second mode dotted in. It was not possible to dot in a probable 
curve for the first mode because no information was available about 
the turning point for this mode. The bimodal data can be represented 
by two terms as in Eq. (18a). Both modes exhibit straight line behavior 
in the large r region, so the RPL (Model 2), was chosen as being suitable 
for both modes. The model fitted takes the form 


n(r) = p 


(r/p 2 ) 


P3" 1 


P 5 <r/P 6 ) 


P 7 “l 


p 3 P 4 


[ p 2 (l+(r/p 2 ) ) 


P 7 P 8 


P 6 (l+(r/p 6 ) ) 


where p 2 , py and p^ are parameters for the first mode, and p,-, p^ , p ? , 

and pg are parameters for the second mode. 

Model 2: p' could not be determined from the data 

3 

for either mode. The DT was overlaid on the p^- CG (Fig. 2A.4) and 

-2 2 

translated Vertically and horizontally until the point (10 , 2 x 10 ) 
coincided with the CG point (0.3, 1.0). The data points for the first 
mode lay along the p^ = 2.0 CP and, in the absence of any information 
about p^ , the values p^ = 3.0, p^ = 2.0 were taken as the best estimates 
of these parameters for the first term. 

To determine pg for the second mode, the DT was 

2 

translated so that (35.0, 3.0) coincided with the CG point (10 , 10). 

The data points, for the second mode lay along the ■ = 2.0 CP, and 
since this CP fitted both the peak and the large r portion of the 
data curve,. p 7 = 3.0, Pg = 2.0 were taken as the estimates. 

To confirm p^ for the second mode, the DT was over- 
laid on the Pg- CG (Fig. 2A.3) and translated vertically and hori- 
zontally until the point (45.0, 2.5) coincided with the CG point 
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(10 , 10 ). The data points in the region of the peak for the second 

mode lay along the p^ = 3.0 CP and confirmed the p^ estimate. 

The next step was to determine p ' and p ' . The DT was 

2 6 

-2 2 

overlaid on the p - CG (Fig. 2A.1) so that the point (10 , 5 x 10 ) 

_2 

coincided with the point (10 ,1.0) on the CG. The data points for 

the first mode lay between the p 2 = 0.1 CP and p 2 = 0.5 CP and in the 

absence of any information about the mode radius, p^ = 0.1 was taken 

as the estimate for the first mode. (The value 0.27, which is the 

smallest r for which data are available, is too high since there is 

no sigh of any curvature at this point.) 

To determine p^ for the second mode, the DT was trans- 

_ 2 

lated vertically until the point (10 , 3.0) coincided with the point 

-2 

(10 , 10) on the CG. In this case, interpolation was necessary since 

the data lay between the p 2 = 1.0 and p 2 = 3.0 CPs. Using Eq. (14b) 
with r^ = 1.6 gave Pg = 3.58 for the second mode. 

The (p^, p^) values for the CG differ from (p^, p^ and (p^, 

Pg) ; thus, the p 2 and p^ estimates had to be adjusted by using Eq. (14c). 
Substituting in this equation gave p 2 = 0.06 for the first mode and 
pg = 2.04 for the second mode. 

The final step was to determine the relative proportions 
of the two modes by determining the scaling factor p^ for each mode. 

Since the peak values of the two modes differ by several orders of 
magnitude, it was decided to determine the scaling factors independ- 
ently, by calculating them at radii where only one mode made a signif- 
icant contribution to n(r) and using Eqs. (18d and 18e) to determine the 
scaling factors a£, a^ for each individual mode. This gave aj = 


1225 for 
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the first mode (at r = 0.27) and = 0.85 for the second mode (at r = 

4.5). Because of the form which the model fitted by the NLLS program 

takes, a j = pj, and a£ = p^ X p^, and thus p^ = 6.9 x 10~\ Thus, 

the initial estimates for this model are p^ = 1225.0, p£ = 0.06, p^ = 

-4 

3.0, p^ = 2.0, p,l = 6.9 x 10 , p^ = 2.04, pj = 3.0 and p^ =2.0. The 
results of the NLLS fit using these estimates are given in Fig. 19. 

CONCLUDING REMARKS 

The uses of the catalog are two-fold. Firstly, it provides a 
catalog of the shapes of the different distributions, illustrating such 
properties as the locations of mode radii, rates of fall-off and poly- 
dispersity. By providing a means for comparing distributions, it aids 
in selecting the model(s) most likely to give a good description of the 
experimental data to be fitted. Secondly, the catalog provides a means 
of estimating the likely values of the model parameters. Often these 
estimates will then be used as initial estimates for a nonlinear least- 
squares or other optimization code. 

The catalog has been used successfully in fitting tropospheric and 
stratospheric aerosol data, bimodal data, . and COSD data. It is importnat 
that the model chosen to fit the data be an appropriate one if useful 
results are to be obtained. In deciding whether a model is appro- 
priate, mode radius, rates of fall-off, and polydispersity should be 
considered. Models which show different behavior for these properties 
from that of the data should be avoided. 
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As the examples showed, there are often several models which will 
adequately describe the data, although some may be better than others. 
Experience with the catalog has indicated that where there is a known 
relationship between a property of the model such as rate of fall-off 
and a particular parameter, that relationship should be used to get 
the best possible estimate for the parameter. 
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SUMMARY TABLE FOR SIZE DISTRIBUTION MODELS 
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Power Law 


Model 2 
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SUMMARY TABLE FOR SIZE DISTRIBUTION MODELS (CONTINUED) 


Model 4 
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Model 6 
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SUMMARY TABLE FOR SIZE DISTRIBUTION MODELS (CONTINUED) 
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Generalized Distribution 
Function 
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Power Law - Generalized 
Distribution Function 
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Parametrized Graphical Catalog 
of Size Distribution Analytic Models 


(For List of Symbols and Acronyms see page vii-viii) 

• Annotation on logarithmic axes gives the power of 10 

. For semi-log graphs ( Figures 6A.1-6C.3, 7A.1-7C.3) 
the annotation on the linear axis (X-axis) is the 
actual value. 

. Blank log-log and semi-log graphs are reproduced on 
pages 156 and 157, respectively. 

. Two transparencies of the blank log-log and semi-log 
graphs are attached at the end of the catalog. 
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Figure 1A.1. Model 1 (power law) for 
n(r). Parameter Set 1.1: P 2 vari- 
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Figure 1B.1. Model 1 for n L (r) 
Parameter Set 1.1. 
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Figure 1C.1. Model 1 for N(>r). 
Parameter Set 1.1. 
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Figure 2A.1. Model 2 (Regularized 
Power Law) for n(r) . Parameter 
Set 2.1: P 2 variable, P 3 = 2.0, 
P 4 = 3.0. 
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Figure 2A.2. Model 2 for n(r). 
Parameter Set 2.2: p 2 variable, 
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Figure 2A.3. Model 2 for n(r). 
Parameter Set 2.3: po variable , 
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Figure 2A.4. Model 2 for n(r). 
Parameter Set 2.4: variable, 
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Figure 2B.4. Model 2 for n L (r) . 
Parameter Set 2.4. 
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Figure 2C.1. Model 2 for N(>r). 
Parameter Set 2.1. 
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Figure 2C.2. Model 2 for N(>r). 
Parameter Set 2.2. 
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Figure 2C.3. Model 2 for N(>r). 
Parameter Set 2.3. 
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Figure 2C.4. Model 2 for N(>r). 
Parameter Set 2.4. 
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Figure 3A.2. Model 3 for n(r). 
Parameter Set 3.2: P3 variable , 
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Figure 3A.3. Model 3 for n(r) . 
Parameter Set 3.3: P3 variable, 
p 2 = 2.0, p 4 = 2.0. 
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Figure 3B.3. Model 3 for n^(r) 
Parameter Set 3.3. 
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Figure 3C.1. Model 3 for N(>r). 
Parameter Set 3.5: P 2 variable, 
P 3 = 20.0, p^ = 0.5. 
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Figure 4A.2. Model 4 for n(r). 
Parameter Set 4.2: P3 variable, 
p 2 = 2.0, p„ = 0.5. 
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Figure 4A. 3. Model 4 for n(r). 
Parameter Set 4.3: p 3 variable , 
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Figure 4A.4. Model 4 for n(r) . 
Parameter Set 4.4: variable, 

P„ = 2.0, p 3 = 5.0. 
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Figure 4B.2. Model 4 for n L (r) . 
Parameter Set 4.2. 
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Figure 4B.4. Model 4 for n L (r). 
Parameter Set 4.4. 
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Figure 4C.1. Model 4 for N(>r). 
Parameter Set 4.1. 
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Figure 4C.2. Model 4 for N(>r) . 
Parameter Set 4.2. 
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Figure 4C.3. Model 4 for N(>r) . 
Parameter Set 4.3. 
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Figure 4C.4. Model 4 for N(>r) . 
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Figure 5A.1. Model 5 (Log Normal 
Distribution) for n(r) . Parameter 
Set 5.1: P 2 variable, p 3 = 0.5. 
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Figure 5A.2. Model 5 for n(r). 
Parameter Set 5.2: P2 variable, 
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Figure 5A.3. Model 5 for n(r). 
Parameter Set 5.3: p 3 variable. 
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Figure 5B.1. Model 5 for n L (r) . 
Parameter Set 5.1, 
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Figure 5B.2. Model 5 for n^(r). 
Parameter Set 5.2. 
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Figure 5B.3. Model 5 for n L (r). 
Parameter Set 5.3. 
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Figure 5C.1. Model 5 for N(>r). 
Parameter Set 5.1. 
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Figure 5C.2. Model 5 for N(>r). 
Parameter Set 5.2. 
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Figure 5C.3. Model 5 for N (>r) . 
Parameter Set 5.3. 
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Figure 6A.1. Model 6 (Normal 

Distribution) for n(r) . Parameter 
Set 6.1: p 2 variable, p^ = 0.5. 
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Figure 6B.2. Model 6 for n L (r). 
Parameter Set 6.2. 


118 



MODELR 

O P2 = 2.50 
u P3 : .300 
PI = 7.70 
RM = 2.50 

n P3 = -500 
PI = 4.68 
RM = 2.50 

/\ P3 = .750 
v PI = 3.19 
RM = 2.50 


PI = 2.46 
RM = 2.50 

K P3 = 1 .50 
PI = 1 .76 
RM = 2.50 

r\ P3 = 2.00 
PI = 1 .42 
RM = 2.50 

q P3 = 3.00 
PI = 1.10 
RM = 2.50 


6 


Figure 6B.3. Model 6 for n^(r). 
Parameter Set 6.3. 
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Figure 6C.1. Model 6 for N(>r). 
Parameter Set 6.1. 
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Figure 6C.3. Model 6 for N(>r). 
Parameter Set 6.3. 
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Figure 7A.1. Model 7 (Generalized 
Distribution Function) for n(r) . 
Parameter Set 7.1: p_ variable. 
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Figure 7A.2. Model 7 for n(r) . 
Parameter Set 7.2.: p variable , 
p 3 = 1.0. 
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Figure 7A.3. Model 7 for n(r) . 
Parameter Set 7.3; p variable, 
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Figure 7B.3. Model 7 for n (r) . 
Parameter Set 7.3. L 



128 



O 


□ 


o 


A 






Q 


Q 

0 


MODELR 

7 

P2 

= 

100. 


P3 

z 

.500 


PI 

• — 

50.5 


RM 

= 

2.30 


P2 

- 

10.0 


PI 

- 

5.50 


RM 

= 

1.15 


P2 


7.50 


PI 


4.25 


RM 

= 

1 .01 


P2 

- 

5.00 


PI 

“ 

3-00 


RM 


.805 


P2 


2.50 


PI 

z 

1 .75 


RM 

z 

.458 


P2 

= 

1 .00 


PI 

z 

1 .00 


RM 

= 

0. 


P2 

- 

o. 


PI 

r 

.500 


RM 

=i 

•000E- 

-02 

P2 

- _ 

.500 


PI 

— 

.250 


RM 

=i 

•000E- 

-02 

P2 


• 950 


PI 

=2 

•500E- 

-02 

RM 

= 1 

•000E- 

-02 


Figure 7C.1. Model 7 for N(>r) . 
Parameter Set 7.1. 
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Figure 1C. 2. Model 7 for N(>r). 
Parameter Set 7.2. 




NC >R 3 



NCR) , 



A P2 
PI 


N P2 
PI 
RM 

A P2 
PI 


Q p 2 
PI 


5.000E-04 

117. 

8.250E-02 

1 .000E-04 
1 .306E+03 
3.750E-02 

5.000E-05 
3 .695E+03 
2.500E-02 

1 .000E-05 
4.111E+04 
1 .250E-02 

: 1 .000E-06 
:5 • 1 39E + 05 
:1 . 000E-02 

1 .000E-08 
9.921E+05 
1 • 000E-02 


-r o i 

RADIUS R > iu.ni 


Figure 8A.1. Model 8 (Power Law — 
Generalized Distribution Function) 
for n(r). Parameter Set 8.1: 

P 2 variable, P 3 = 40.0, p^ = 2.0. 
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Figure 8A.2. Model 8 for n(r) . 
Parameter Set 8.2: p variable, 
-4 J 

d. = 5.0 x 10 , p. =2.0. 
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Figure 8A.3. Model 8 for n(r) . 
Parameter Set 8.3: p. variable, 
P 2 = 5.0 x 10 -4 , p_ = 40.0. 
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Figure 8B.1. Model 8 for n (r) . 
Parameter Set 8.1. L 
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Figure 8B.2. Model 8 for n^(r) . 
Parameter Set 8.2. 
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Figure 8B.3. Model 8 for n (r) . 
Parameter Set 8.3. 
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Figure 8C.1. Model 8 for N (>r) . 
Parameter Set 8.1. 
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Figure 8C.2. Model 8 for N(>r). 
Parameter Set 8.2. 
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Figure 8C.3. Model 8 for N(>r) . 
Parameter Set 8.3. 
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Figure 8D.1. Model 8B (Power Law-- 
Generalized Distribution Function) 
for n (r) . Parameter Set 8.4: 
p 2 variable, p 3 =40.0, p^ = 4.0. 
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Figure 8D.2. Model 8B for n(r). 
Parameter Set 8.5: P3 variable, 
P- = 5.0 x 10 -4 . o. = 4.0. 


NtR) 


$ 

2 p 

o -2 


-3 


-4 




MODELR 8B 

P2 =5 -000E-04 
P3 = 40.0 
P4 = 1 .50 
PI = 5.01 
RM = .200 

P4 = 2.00 
PI. = 12-6 
RM = .150 

P4 = 2.50 
PI = 35.7 
RM =1 .000E-01 

P4 = 3-00 
PI = 117. 

RM =8 .250E-02 

P4 = 3.50 
PI = 451 . 

RM =5 .500E-02 


P4 = 4.00 
PI =2 -287E+03 
RM =2 .750E-02 

P4 = 4.50 
PI =1 .524E+04 
RM =2 • 000E-02 

P4 = 5.00 
PI =1*1 36E+05 
RM =1 .750E-02 

P4 = 6.00 
PI =7 .434E+06 
RM =1 -500E-02 

P4 = 7.00 
PI =5 • 794E+08 
RM =1 .250E-02 


~>l i I., I, I. l-l 111 J i. ,1 — 

-2 -1 0 


2 


RRDIUS R, jum 


Figure 8D.3. Model 8B for n(r). 
Parameter Set 8.6: p 4 variable, 
p 2 = 5.0 x 10 -4 , p, =40.0. 
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Figure 8E.1. Model 8B for n L (r). 
Parameter Set 8.4. 
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Figure 8E.2. Model 8B for n L (r) • 
Parameter Set 8.5. 
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Figure 8E.3. Model 8B for n L (r) . 
Parameter Set 8.6. 
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Figure 8F.1. Model 8B for N(>r) . 


Parameter Set 8.4. 
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Figure 8F.2. Model 8B for N(>r). 
Parameter Set 8.5. 
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Figure 8F.3. Model 8B for N(>r) . 
Parameter Set 8.6. 
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Figure 11. Example of fit to Junge 
data using Model 2. (Regularized 
Power Law) . 
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Figure 12. Example of fit to Junge 
data using Model 4. (Inverse 
Modified Gamma Distribution) . 
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Figure 13. Example of fit to Junge 
data using Model 5. (Log Normal 
Distribution) 
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Figure 14. Example of fit to Junge 
data using Model 8B. (Power Law- 
Generalized Distribution Function) 
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Figure 16. Example of fit to Alaska 
data using Model 3. (Modified 
Gamma Distribution) . 
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